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Abstract. A simple mechanical system consisting of a mass and a dry friction element is 
considered. The mathematical model of such mechanical system is described. We show that 
the considered system corresponds to a skew product over a piecewise smooth mapping of a 
segment. We demonstrate that this base map for some value of parameters a robust chaotic 
dynamics can be observed. 
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Introduction. 

Systems with dry friction form a wide class inside discontinuous dynamical systems. They 
appear in many applications, especially in manufacturing systems: vibrating conveyors, per- 
cussion drilling, metal cutting, etc (see [2|3Pf5|6|12fl3fl4fl6fl9f22f23 ] and references therein 
for review). Their properties manifest many principle differences with ones of smooth dy- 
namics. For instance, the uniqueness theorem is not valid any more. An approach to study 
such systems has been developed by A. F. Filippov [7]. He offered to consider piecewise 
continuous systems of differential equations as families of vector fields, defined on disjoint 
domains of the phase space and define auxiliary tangent flows on boundaries, respecting 
limit directions of vector fields. This approach reduces a discontinuous system to differen- 
tial inclusions. Moreover, the phase space may become multidimensional i.e. a set of initial 
data of the full dimension may be transferred to a set of a lower dimension. The theory of 
discontinuous systems and specific bifurcation is well-developed [2|3f5|7p0fl4lll6lll8f22j . It 
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is also well-known that chaotic dynamics frequently occurs in such systems, particularly in 
ones with dry friction [T]IMK]l6]IMM71[T9]l20^ 




Apart from numerical and experimental simulations, the most common analytic approach 
involves a reduction of dimension. For some systems with a dry friction it is possible to 
demonstrate that there exists an invariant set of dimension 1 where the attractor resides. 
A method to find this attractor has been developed by M. Wiercigroch, E. Pavlovskaya and 
A.Krivtsov in papers [T2|T3] and, in its general form, in the paper [iTj. In our paper we use 
some ideas of this approach. Another powerful method has been proposed by R. Szalai and 
H. M. Osinga [20j. They have proved that for a general class of systems with a dry friction 
that the attractor resides in a polygon type set and demonstrated a possibility of a chaotic 
dynamics there. Later [21] they have shown, using a modification of their method that some 
complex structures like Arnold tongues can be observed in a neighborhood of the so-called 
grazing-sliding bifurcation [2|3p8] . 




The main aim of this paper is to provide a new method which allows to find chaotic invariant 
sets in systems with a dry friction. To demonstrate this method, we use a very simple 
example of a system with dry friction, first considered be A. Krivtsov and M. Wiercigroch 
fn\ . First of all, we show that the considered dynamical system engenders a discontinuous 
mapping of a segment. Here we use ideas from [17]. Then we study properties of this mapping 
which allows us to find two disjoint segments such that the image of every one of them covers 
their union. Moreover, we prove that the considered mapping is continuous on the union of 
these segments. This allows us to apply well-known techniques of one dimensional dynamics 
[TT| Part 3, Section 15] and to demonstrate that a kind of chaotic dynamics, similar to one 
described by T. Li and J. Yorke [15] is there. 

The main advantages of the offered method are the following. 

(1) We can obtain chaotic sets which, in general, are not attractors. 

(2) For simple systems with dry fiction, the offered method gives coefficient type criteria 
of chaos. 

(3) Though we need a presence of a small parameter in our proofs, it is possible to estimate 
numerically how small this parameter must be. The presence of chaotic invariant sets 
does not correspond to a neighborhood of any bifurcation. 

(4) The presence of chaotic behavior (but not corresponding invariant sets) is robust. 

(5) A corresponding invariant measure can be described using techniques of [HI Part 3, 



The paper is organized as follows. In Section 1 we introduce the mathematical model of 
the considered system and describe possible regimes of motion. In Section 2 we define the 
main object of our investigation: the 1 dimensional mapping, corresponding to phases of 
switching for solutions. For this, we describe all possible scenarios of behavior of solutions. 
In the next section we describe the main properties of the introduced mapping. The main 
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part of the proof is given in Section 4 where we study how the segments of continuity of 
the constructed mapping look hke and find out two segments of continuity whose images 
cover their union. Section 5 is technical. We would like to prove existence of periodic points 
of all possible periods but we cannot apply the theory of mappings of a segment directly 
since we deal with a discontinuous mapping. However, still we can use standard methods 
of this theory to finish our proof. The main result is formulated in Section 5. A discussion, 
including notes on the robustness of the obtained set and some plans on the future research 
is given at Section 7. 

1. Description of the mathematical model. 

Consider a single degree-of-freedom mechanical system, consisting of a point mass and a 
delimiter with dry friction (Fig. 1) which gives a simple model of percussion drilling. This 
system consists of a unit mass, whose motion is controlled by a harmonic external force F{t) 
which is a sum of a positive constant component equal to 2b and a harmonic component 
of a positive amplitude a and a period equal to 27r. Also, the considered system includes 
a delimiter which provides an additional dry friction as soon as the mass reaches it. The 
maximal value of this friction is q. Here and later we always suppose that all considered 
parameters are non-dimensional. 

Our main aim is to prove that provided some additional conditions are satisfied, the dynam- 
ics of the considered system is robustly chaotic in topological sense. Let x be the current 
position of the mass and y be one of the delimiter. We assume that the inequality x < y is 
always satisfied i.e. the mass cannot penetrate through the delimiter and that the delimiter 
cannot move to the left so the value y is always non-decreasing. 




Fig. 1. The considered mechanical system. 



Consider the value 'do e [0,7r/2] such that 



71-^0 = cot(V2)- 




This value is unique and t^q fa 0.81047, sini?o 



0.724611. 



We make the following assumptions on the parameters of the system: 



a > 0, 



6> 0, 



q e (asini?o, a). 



(2) 
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We shall always suppose that b <^ a (which implies b <^ q). This means that we always 
suppose that the ratio b/a is as small as necessary. 

There are five types of motions of the considered system. 

(1) Free motion (f). This motion takes place if x < y i.e. the mass and the delimiter do 
not interact. Then 

X = = asint + 26; y = 0. (3) 

(2) Progression (p). In this regime we have an additional friction. The motion is defined 
by equations 

X = F(t) — q = asint + 2c = asint + 2b — q; y = x. (4) 

(3) Stop (s). Here the mass and the delimiter are both immobile, i.e. 

X = y, X = y = 0. (5) 

(4) Instantaneous stop (is) . This happens if Condition (|5]) is satisfied for a fixed instant 
of time but is not true in its vicinity. So, this happens if the system switches from or 
to free motion or from/to motion with progression. 

(5) Transition from the free regime to motion with progression (fp). This is an 
instant such that the mass was moving in a free regime before and switches to motion 
with progression there. 

We always suppose that the free motion, motion with progression and the stop regime are 
observed at open intervals of time. This allows us to classify all instants of transition. 

Let to < ti be zeros of the function F{t) — q and t2 < ^3 be ones of the function F{t) on 
[0, 27i]. Later on we consider the phase ip = t mod 2tt. Here ip E = M/27rZ. 

Now we describe how we may proceed from one regime to another. 

If the mass, moving according to Eq. collides with the delimiter with a non-zero velocity, 
it proceeds to the progression regime. The transition (f) — t- (is) — )■ (s) is possible at t = 6'o 
if and only if x{9o) = y{0o), x{9o) = 0. If F{9o) > then x{t) is negative before t = 9o which 
is impossible. If F{6o) < 0, the mass returns back, "ignoring" presence of the delimiter. We 
have an instantaneous stop there. Otherwise, F{6q) = 0. If F{6o) < then again we have 
x{t) < in a left neighborhood of Oq. So, F{Oq) > and, consequently, = h mod 2tt. 
Then we have an instantaneous stop and the stop regime later on. 

The progression regime is characterized by condition i > 0. As soon as the solution x{t) 
reaches its maximal value, corresponding to t = 6', the mass stops and then proceeds to 
free motion if F{t) < in a right neighborhood of 6'. Also, we may proceed to the "long" 
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(non-instantaneous) stop regime if F{t) G [0, q] there. Otherwise, we stay in progression 
regime. Anyway an instantaneous stop takes place. 



The system stays in stop regime until we have F{t) < or F{t) > q. In the first case we 
proceed to the free regime, in the second one we start progression. See Fig. 2 for illustration. 




Fig. 2. Possible regimes of the system and transitions according to the phase. 
2. Reduction to dimension 1. 

Solutions of Eqs. ([3]), @ and ([5]) can easily be written down. If x{9o) = xq, x{9o) = xi we 
have 

x{t) = —a sin t + b{t — OqY + {xi + a cos 6*0) (t — 60) + xq + a sin Oq (6) 
for Eq. (E]) (free motion) and 

x{t) = —a sin t + c(t — OqY + {xi + a cos 6*0) (t — 9q) + xq + a sin Oq (7) 

for progression. For stop regime we always have xi = and x{t) = xq. 

Let us show that none of the considered regimes can be eternal. 

Starting from the free motion, the mass will always return to the delimiter since 6 > 0. Then 
we immediately proceed to progression motion unless we reach the delimiter with a velocity 
equal to zero. Let this zero velocity approach happens at the instant t which corresponds to 
the phase r. Note that r G (^2, ^3] otherwise the velocity would be negative in a left vicinity 
of r. Then if r G (^2,^3) the mass stops and starts moving in the free regime immediately 
after that. If r = ts the mass stops until the next instant to + 2nk. 

In the motion with progression the derivative x vanishes soon or later since c < 0. If this 
happens when (p G [^2,^3) we immediately proceed to free motion. Otherwise, the mass 
stops. If this happens for ip G [^1,^2) the mass stops until if = t2 and then switches to free 
motion. Note, that progression cannot be stopped while t G [to, ti). If it is stopped on [^3, to) 
the mass waits the next instant to + SvrA; and then starts moving according to Eq. (|4]). In 
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this case, we have 6^0 = to and xi = in Eq. ([7]). Consequently, 

x{t) = — a(cost — cos to) + 2c(t — to)- (8) 

This function increases (and, therefore, cannot vanish) until t = ti. However, since b <^ q, 
we may say that costs > cos to and the right hand side of Eq. ([8]) is negative for t = t^. So, 
the motion stops somewhere at [ti,t3) and then proceeds to the free flight regime. The stop 
regime may be finished by a transition to a free motion at t = t2 + 27ck or by transition 
to the motion with a progression at t = to + 2'n'k. It is impossible neither on (to,ti) nor on 

(t2,t3). 

All together, this means that wherever and whenever the motion starts, finally, it has a 
transition to a free motion via an instantaneous stop. 

Take an initial instant 6 of such transition. There we have 6 G [t2,t3], x{6) = y{6) (without 
loss of generality, we may assume that this value is zero) and x{6) = 0. Then the value 6 
uniquely defines the farther dynamics. 

If T{6) is the phase of the next transition to free motion from stop or progression, it is 
uniquely defined by 6. So, we may consider the ID mapping T : [t2, t^) O which is, generally 
speaking, discontinuous. Considering this mapping, we lose some information about initial 
dynamical system, for instance, we do not know any more how the delimiter is shifted. 

One of the following scenarios takes place for the considered motion. 

(1) Scenario A: (f) ^ (fp) ^ (p) ^ (is) ^ (f); 

(2) Scenario B: (f) ^ (fp) ^ (p) ^ (is) ^ (s) ^ (is) ^ (f); 

(3) Scenario C: 

(f) ^ (fp) ^ (p) ^ (is) ^ (s) ^ (is) ^ (p) ^ (is) (s) ^ (is) ^ (f); 

(4) Scenario D: 

(f) ^ (fp) ^ (p) ^ (is) ^ (s) ^ (is) ^ (p) ^ (is) ^ (f). 

Also, there are two degenerate scenarios, corresponding to a zero-velocity stop of the free 
motion corresponding to t = ts + 2'n'k: 

(1) Scenario C: 

(f) -> (is) ^ (s) ^ (is) ^ (p) ^ (is) ^ (s) ^ (is) ^ (f); 

(2) Scenario D': 

(f) ^ (is) ^ (s) ^ (is) ^ (p) ^ (is) ^ (f); 
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Here we omit all possible instantaneous stops after which the motion returns to the same 
regime. This does not hurt to equations of motion. However, such stops play an important 
role since they correspond to discontinuities of stroboscopic mappings. Let us study them 
more carefully. 



3. Properties of the mapping T. 



Let us introduce two auxiliary mappings Ti and T2. The first one corresponds to the first 
after 9 approach to the delimiter. Denote the corresponding instant by 9i, i.e. Ti{9) = 9i. 
Note that the image of Ti is an instant, not phase, so it can be greater than 27r. 

The value 9i corresponds to the first zero of the equation 

b{t -9y/a = smt- cos 9{t - 9) - sin 9. (9) 
satisfying the condition t > 9. 

The left hand side of Eq. is always positive and proportional to the small parameter 
b/a. The right hand side is initially positive and grows faster than the left hand side. This 
is the distance between the graph of sine function and the tangent line to it, drawn at 9. 
Note that sin^ < —b/a < 0. If ^ G (37r/2,27r) i.e. cos^ > 0, the graph and its tangent line 
intersect once again on {9, +00) and, therefore 9i — 9< 27i. Otherwise, they do not intersect 
and, for small b, there exists a constant C > which does not depend on 9 and is such that 

^i-^> -=^cos^. (10) 


The mapping Ti is, in general, discontinuous. All possible discontinuities correspond to the 
case when 9i is not a simple zero of (Q. Namely, derivatives of the left and the right hand 
side of coincide for t = 9i which means that the following condition is satisfied 

— a cos 9i + 2b{9i — 9) + a cos 

^^^/sin^i — ^ \ ^ Q ^^^^ 



4. Points of discontinuity. 

Lemma 1. The intersection of the set of discontinuity points of the mapping Ti with Lq = 
[10l7r/100,t3] ^s finite. 

Proof. If b/a is small ta < lOlvr/lOO. Note that if 

9^ < 9^, 9^'^ e Lo, (12) 
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then Ti{6^) > Ti{9'^). Indeed, consider the function 

hit, 6) = b{t - ef/a - sin t + cos e{t - 6) + sin 6. 
The derivative of this function with respect to 9 equals to 

- {^ + sm9^ (t-e) 

and is positive if ^ G Lq, b/a is small and t > 6. li 9^ and 9'^ satisfy (fT2|) and is such 
that 9^) = then h{t^, 9"^) > and the function h{t, 9"^), negative in a right vicinity of 
1 = 9"^ must have a zero on {9^,t^). 

So the function Ti is monotonous. Note that if 9^ G Lq is a point of discontinuity and 
ti = Ti(^i) then 

(fh, , 2b 



-{t\9^) = — + sint^ 



Note that due to ( |T0|) and (ITTl) we have sint^ = sin6'^ + 0{b/a) which implies that this 
second derivative must be greater that a fixed negative value. So there exists a p > such 
that for all discontinuity points corresponding "jumps" of the function Ti must be greater 
than p. This proves that the number of discontinuity points on Lq is finite. □ 

It follows from Eq.llTTl) that for any discontinuity points 9 of the mapping Ti we must have 

cos9i = -cos9 + 0{{9i-9)-^). (13) 

Another useful note is that Ti{9) mod 27r G [^2,^3) cause we cannot have a transition to a 
free regime in other points. Meanwhile, for these discontinuity points, the second derivative of 
solutions (right hand side of Eq. ([3])) must be negative, since x < y (delimiter is always to the 
right of the mass). So, the corresponding value 9i must belong to [37r/2,t3] if 6* G [t2,37r/2]. 

In non-degenerate scenarios ((A)-(D)) at the moment t = 9i progression regime starts. The 
initial velocity of the motion is xi = —a cos 9i + 2b{9i — 9) + a cos 9. The dynamics of this 
velocity is described by the formula x{t) = Xi + 2c{t — 9i)~acost + acos{9i). The progression 
regime stops as soon as this derivative becomes negative and the next transition ^2 = ^2(6') 
to free flight or to the stop may be found from equations 

2b{9i -9)+acos9 + 2c(^2 - 9i) - a cos ^2 = 0. 



Lemma 2. There exists a bo{q) > such that if b < bo{q) the following statement is true. 
Starting progression at the point to with an initial velocity equal to zero, the motion must 
stop at the instant t^ G [^1,^2]- Consequently, Scenarios D and D' are impossible. 
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Proof. If x{t) is a solution of Eq.Q with x{to) = 0, we have 



X{t2) 



a cos t2 + a cos to + 2c(t2 — ^o) • 



(14) 



In order to prove that this x{t) vanishes somewhere at [to ^2) it suffices to prove that the 
right hand side of Eq. f|T^ is negative. Instead of this one could prove that 



Here we replaced t2 with vr in (HM and respected the fact that 2c = 2b — q. If we demonstrate 
for a fixed q that 



then for this fixed q there exists a 60 (O') > such that if 6 < &o(9) then the estimate (fT4|) is 
true. 

Inequality (|T5i) is equivalent to the following: cot (to/2) < tt — to which is true if to > ^0 (see 
Eq.([T])) or, equivalently if g > sini^oct. □ 

Lemma 3. The map T is such that T{6) = 62 if O2 G [t2,t3). Otherwise, if a motion with 
stop has been observed (Scenarios B and C), T{9) = t2. 

Proof. If 62 G [t2,t3) then a motion, starting near the delimiter with the velocity, equal 
to zero, corresponds to the free regime, so T = 62- If O2 G [t3,to], the mass stops until the 
instant to then starts moving in progression regime until t = t4, stops until t2 and proceeds 
to a free regime. If 62 G [ti,t2) the motion stops until t = t2 and also proceeds to the free 
regime. Since 62 cannot belong to (to,ti), we obtain the statement of lemma. □ 

To finish our proof, we need the following lemma. 

Lemma 4. There exist two disjoint subsegments Jo o-nd Ji of the segment [t2,37r/2] such 
that T{Ji) D [t2,3n/2] and T is continuous on both segments Ji. 

Remark. Then we may say without loss of generality that T(Jj) = [t2, 37r/2]. 

Proof. Let Li be the arc [1977r/100, ts] and L2 be the arc [IOItt/IOO, Slvr/SO]. The first arc 
is correctly defined if the ratio b/a is sufficiently small. 

Let 9 ^ Li. Consider the solution x{t) such that x{9) = 0, x{9) = 0. Then, direct calculations 
show that for 6 = 0, the considered solution starts from a short period free motion and then 
stops after a short period of motion with a progression. So, i:(to) = 0. Due to continuous 
dependence of the solution on the parameter b, the same is true provided the ratio b/a is 
sufficiently small. In this case, as we have already proved T{9) = t2. 



a + a cos to — (g — 26) (tt — to) < 0. 



a + a cos to — g(vr — to) < a(l + cos to — sinto(vr 



to)) < 0. 



(15) 



9 



Let zi < Z2 < ■■■ < Zn {n > 0) be discontinuity points of the map Ti inside the interval L2. 
Denote Zo = IOItt/IOO, = Slvr/SO. 

Suppose that n < 3. Then there exists i G {0, . . . ,n} such that there is a subsegment 
X C {zi, Zj+i) of the length not less than 7r/500. Let 9 E Z, 9i = Ti{6). 



It follows from Eq. ([9]) that 



a ( ^ sin ^1 — sin 6'' 

"=6 e,-e 



Consequently, if h/a is sufficiently small, values {9i mod 2t{ : 9 eX} cover [0, 27r). Since Ti 
is continuous on X, we have 

2h 

— (^01-9) ^ cos 01 - cos ^ 
a 

or, equivalently, 

sin 6'i — sin 6^ , ^ ^, , . , 

2 -i (cos 01 + cos 0) ^0 (16) 

9i — 9 

everywhere on X. However, due to Eq. ( |T0|) . the maximum of the left hand side of in- 
equality ( !T6|) is positive while the minimum is negative. So this inequality cannot hold true 
everywhere. 

So, n > 3. Then it suffices to prove that T{{zi, Zi+i]) D [t2,37r/2] for alH = 1, . . . , — 1. 
Note that estimate f lT3|) implies that Ti(zj_|_i) > 37r/2 if b/a is small. Since x{Ti{zi+i)) = 0, 
the corresponding motion proceeds to the free flight immediately after t = Ti (2:4+1) and, 
consequently, T{zi+i) = Ti(2j+i) G Li. 

On the other hand, Ti{9) — y Ti{Ti{zi)) as 9 — > Zi + (in the limit case, we have a motion, 
which "ignores" presence of the delimiter). Since Ti{zi) G Li the corresponding motion is 
in the stop regime for t = to and, consequently, T{9) = t2 for all 9 from the right vicinity of 
Zi. This finishes the proof. □ 



5. Infinite set of periodic points. 



So, for the mapping T we have obtained two disjoint segments Jq and Ji which are subsets 
of the arc [^2, ^3] of the unit sphere such that for both i = 0,1 mappings T\j. are continuous 
and T{Ji) D JqU Ji- Then we prove that for any m G N the mapping T has a point of the 
minimal period m. The prove is absolutely standard, but it has to be given, since formally, 
the existing theorem is applicable to continuous maps of a segment. Also, note that we do 
not claim that periodic point of a fixed period is unique. There might exist a segment (or, 
even a Cantor set), consisting of such points. Take a sequence {cjfc G {0, 1} : G NU{0}}. 
First of all, we prove that there exists a point p G J^q such that T^{p) G J^^. for any /c G N. 
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There exists a segment Jo-oo-i C Ja^ such that T( Jo-go- J = Jo^. Then, we may find a segment 
Joocricr2 <^ •^o-ocri sucfi that T(Jo(,cti(72) = ^02- Repeating this procedure, we obtain a nested 
sequence of segments 

The corresponding intersection is non-empty and, consequently contains a desired point p. 

Fix a number m and consider the sequence cr, obtained by infinite repetition of a finite 
sequence 0, . . . , 0, 1 of the length m. Let / = J0...01 (see above). Then / C T"^{I) and the 
mapping T"^ is continuous on this segment. Applying Weierstrass principle to the continuous 
function T™(x) — x on the segment I, we obtain a periodic point. Clearly, this point cannot 
be one of a lower period. 

6. Conclusion. 

Let us formulate the principle result of the paper as a single theorem. 

Theorem 1. For all a and q, satisfying inequalities ([2]) there exists a bo = bQ{a,q) > 
such that for all b G (0, bo) the mechanical system, described by equations ([3]), @ and ([5]) is 
chaotic in the following sense. The phase of transition to free flight of the motion uniquely 
defines the phase of the next transition. This defines a discontinuous mapping T from the 
segment [^2)^3] ^^^o itself. Here t2 and t^ are zeros of the function 

F{t) = asint + 6. 

There exist two disjoint segments Jq and Ji of the segment [t2,37r/2] such that T{Ji) D 
[t2,37r/2] andT is continuous on both segments Ji. Particularly, there exists an infinite set 
P of periodic points of the mapping T. Minimal periods of points of P are unbounded. 

7. Discussion and future work. 

First of all, let us note that the obtained chaotic dynamics is robust. Of course, we cannot 
say anything about stability of points of the set P. Every particular point of this set may 
appear or disappear if we slightly change parameters a, b and q. The power of the set P may 
be continuum for some values of parameters, while this set is countable for other values. 
Neither, the method we offer does not specify the topological structure of the set P and one 
of its closure. 

However, we can select a family of segments Jq and Ji from the statement of Theorem 1 so 
that boundary points of these segments locally continuously depend on parameters a, b and 
q. This means that the fact of presence of the one dimensional turbulence in the considered 
system is robust. The same is true for the fact of existence of an infinite set P. Moreover, 
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for all fixed values a, b and q, satisfying inequalities ([2]), there exists an e > such that for 
any function 



such that 



G{t,x,x) : X 



\G{t, x,x)\ < e, 



— ^ 



dG 



< e 



d{t, X, x) 

for all t, X and x an analog of Theorem 1 is true for the system, where equation 
replaced with 

X = F{t) + G{t,x,x); y = 0, 

equation (jlj) is replaced with 



IS 



X 



F{t) + G{t,x,x) - q; y = 0. 



and equation is the same. 



Particularly, the presence of the considered chaotic dynamics must be observed in simula- 
tions and experiments. However, in this paper, we are not going to study the general case. 
We just offer a method how a non-classical chaos may be found. We plan to use this methods 
for more general systems with a dry friction (see as an example) and provide for these 
"real life systems" theoretical results accompanied with simulations and experimental data. 
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